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Abstract 



We analyze the effect of the Sun's gravitational field on a flow of cold 
dark matter (CDM) through the solar system in the limit where the velocity 
dispersion of the flow vanishes. The exact density and velocity distributions 
are derived in the case where the Sun is a point mass. The results are ex- 
tended to the more realistic case where the Sun has a finite size spherically 
symmetric mass distribution. We find that regions of infinite density, called 
caustics, appear. One such region is a line caustic on the axis of symme- 
try, downstream from the Sun, where the flow trajectories cross. Another 
is a cone-shaped caustic surface near the trajectories of maximum scattering 
angle. The trajectories forming the conical caustic pass through the Sun's 
interior and probe the solar mass distribution, raising the possibility that the 
solar mass distribution may some day be measured by a dark matter detec- 
tor on Earth. We generalize our results to the case of flows with continuous 
velocity distributions, such as that predicted by the isothermal model of the 
Milky Way halo. 



PACS number: 98.80 Cq 
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I. INTRODUCTION 



There are compelling reasons to believe that the dark matter of the universe is consti- 
tuted in large part of non-baryonic collisionless particles with very small primordial velocity 
dispersion [|TJ. Generically, such particles are called cold dark matter (CDM). The leading 
CDM candidates are the axion and the lightest supersymmetric partner in supersymmetric 
extensions of the Standard Model of elementary particles. The latter candidate is a weakly 
interacting massive (mass in the 10 to 100 GeV range) particle, or WIMP for short. The 
former has mass in the 10~ 5 to 10 -6 eV range and is extremely weakly interacting, so weakly 
in fact that it was once thought "invisible" . 

Wherever a galaxy forms, the surrounding CDM falls into the galactic gravitational 
potential. Unlike baryons, CDM is collisionless and sloshes back and forth in the potential 
well, producing a galactic halo of size much larger than the extent of the visible matter. The 
local (i.e., on Earth) dark matter density of the Milky Way halo is of order 10 -24 gr/cm 3 . 
This local dark matter is in principle detectable and considerable effort is being spent on 
projects devoted to this goal. Dark matter axions can be detected by their conversion to 
microwave photons in an electromagnetic cavity permeated by a strong static magnetic field 
0. WIMPs can be detected by observing the recoil energy of nuclei off which WIMPs 
have scattered elastically ||£|]. WIMP dark matter may also be detected indirectly by 
observing, on Earth, neutrinos produced in the annihilation of WIMPs which accumulated 
in the interior of the Sun ||. 

From the point of view of these dark matter searches it is desirable to know as much as 
possible about the phase-space distribution f(r, v) of CDM in galactic halos in general, and 
particularly in the Milky Way halo. One wants to know not only the local density but also 
the local velocity distribution. In the cavity dark matter axion detector, the kinetic energy 
spectrum of dark matter axions can be measured with great resolution ||. Knowledge of 
the velocity distribution may therefore be exploited to increase the sensitivity of the search. 
In WIMP searches, knowledge of the velocity spectrum is necessary to calculate the annual 
modulation of the signal . This annual modulation is an important means by which the 
signal to background ratio can be enhanced. 

Yet, in spite of this urgent need of the ongoing cold dark matter searches, there is at 
present no consensus about the phase-space structure of galactic halos. Let us discuss briefly 
the three main approaches which have been put forth. 

The oldest description of galactic halos is the isothermal model f?J. It postulates that the 
phase-distribution is given by Boltzmann's law, the energy of each particle being its kinetic 
energy plus its gravitational potential energy in the presence of all the other particles. It 
can be shown that the dark matter density then falls off as d(r G ) cx 4- for large values of the 
galactocentric distance tq- Therefore the model predicts that galactic rotation curves are 
flat {v (r G ) oc constant) at large r G , and that is indeed what is observed. Another virtue of 
the model is that it predicts that at small r G , within some radius called the core radius a, the 
density d goes to a constant, and hence that the halo contribution to galactic rotation curves 
VhaXoirc) oc tq at small r G . That also is consistent with observation. The model predicts 
all properties of galactic halos in terms of just two parameters: the velocity dispersion y/W 
and the core radius a. In fits to observed properties of our own halo, v^/ 2 is found to be 
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approximately 270 km/s and a of order several kpc. 

The isothermal model is simple and gives an excellent overall description of what is 
observed, but there is little justification for its basic assumption, that the dark matter 
particles in galactic halos are thermalized. Thermalization may have occurred during an 
early phase of 'violent relaxation' || at the time of collapse of the protogalaxy. It is possible 
that the dark matter that fell onto the Galaxy at early times was thermalized by that process. 
But an isolated galaxy such as the Milky Way is still accreting dark matter today because 
it is surrounded by a sea of dark matter. (The presence of M31 and the other members of 
the local group is not relevant for the purpose of this discussion.) There is no evidence of 
violent relaxation occurring in the Galaxy at the present time. The infalling dark matter 
particles form flows which slosh back and forth in the gravitational potential of the galaxy. 
The flows do not get thermalized over the age of the universe ||. Gravitational scattering by 
inhomogeneities, such as globular clusters and molecular clouds, is too weak to diffuse the 
flows entirely, except for those flows which always remained within the inner parts [tq ~ 20 
kpc, say) of the Galaxy. 

The flows can be described in terms of the evolution of the 3-dim. sheet on which the 
dark matter particles lie in 6-dim. phase-space || . Before the onset of galaxy formation, the 
equation of the sheet is v = H(t)r + Av(f,t) where H(t) is the Hubble expansion rate and 
Av(f, t) is the peculiar velocity field associated with density perturbations. The thickness 
of the sheet is the primordial velocity dispersion 5v. The primordial velocity dispersion of 
the leading cold dark matter candidates is extremely small, of order 10~ 12 c for WIMPs and 
3 • 10 _17 c (at most) for axions. The sheet cannot tear and hence its evolution is constrained 
by topology. Where a galaxy forms, the sheet wraps up in phase-space, turning clockwise 
in any two dimensional cut (x, x) of that space, x is the physical space coordinate in an 
arbitrary direction and x is the associated velocity. The outcome of this process is a discrete 
set of flows at any point in a galactic halo. Two flows are associated with particles falling in 
and out of the galaxy for the first time (n = 1), two other flows are associated with particles 
falling in and out of the galaxy for the second time (n = 2), and so on. The number of such 
flows at the Earth's location in the Milky Way is of order 100 ||. The flows of particles 
which have fallen in and out of the Galaxy only a small number of times in the past (n ~ 20) 
are not expected to be thermalized. The flows with larger n may be thermalized because 
they are composed of particles which fell in a long time ago. 

Caustics appear wherever the projection of the phase-space sheet onto physical space 
has a fold [pT0|— 112[] . Caustics are generically surfaces in physical space. On one side of the 



surface there are two more flows than on the other. At the surface, the dark matter density 
is very large. It diverges there in the limit of zero velocity dispersion. There are two types 
of caustics in the halos of galaxies, inner and outer. The outer caustics are topological 
spheres surrounding the galaxy. They are located near where a given outflow reaches its 
furthest distance from the galactic center before falling back in, and are described by A2 
(or 'fold') catastrophes. The inner caustics are rings ||10|| . They are located near where the 
particles with the most angular momentum in a given inflow reach their distance of closest 
approach to the galactic center before going back out. A caustic ring is a closed tube whose 



cross-section is a D_4 (also called elliptic umbilic) catastrophe [12|. The existence of these 



caustics and their topological properties are independent of any assumptions of symmetry. 
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A second approach to describing the structure of galactic halos is in terms of these flows 
and their associated caustics [P, [l^[i^|l4i|l^[i^[T6| . Let us call that approach the 'caustic 
ring model' since caustic rings are one of its prominent features. In the model the caustic 
ring radii, and indeed the whole phase-space distribution, are calculated assuming that halo 



formation is a self-similar process |IT| , |T3| . The model depends on only three parameters: 
the galactic rotation velocity t> ro t, a parameter j max which is proportional to the amount of 
angular momentum in the halo, and a parameter e which describes the profile of the initial 
overdensity from which the galaxy forms, e is related to the slope of the evolved power 
spectrum of primordial density perturbations and hence, to the extent that that spectrum 
is predicted by CDM cosmology, e is not a free parameter. If the primordial spectrum is 
Harrison-Zel'dovich, e is in the ranee 0.2 to 0.35. m. In that range, the self-similar model 



predicts that the density dire) cc T at large r G . Hence the model predicts that galactic 
rotation curves are flat at large r^, which is consistent with observation. When j max = 0, 
the d(r G ) oc 4- behavior extends to small r G - The 4- singularity at r G = is due to the 

r G r G 

fact that the dark matter particles move on radial orbits in that limit. However, the dark 
matter particles are expected to carry angular momentum. When j max ^ 0, the density is 
depleted within an effective core radius r c . As a result, the halo contribution t>haio( r G) to 
the rotation curve goes to zero as r G — > 0, which is also consistent with observation. 

Evidence for the presence of caustic rings in our galaxy was found in the distribution of 
rises in the Milky Way rotation curve, and in a triangular feature in an IRAS map of the 
galactic disk ||16|| . The feature is reminiscent of the shape of the _D_4 catastrophe. When 
the model is fitted to these observations, its parameters are determined to be: v TOt ~ 220 
km/s, and j max — 0.26 for e = 0.3. The model predicts the dark matter flows and caustics in 
the solar neighborhood. Because of our proximity to the fifth caustic ring, one of the flows 
dominates over all others combined. Its density and velocity vector (known only up to a 
two- fold ambiguity) relative to the non-rotating restframe of the Galaxy are |16| : 

4 = 1.7-10- 24 -^, t7 5 = (47O0 G ±lOOf G )— . (1.1) 
cur s 

tq is the local unit vector in the direction away from the galactic center, and 4>g is in the 



direction of galactic rotation. A list of the other flows is given in ref. [fL4j]. The caustic ring 
model describes the minimum phase-space structure that a galactic halo must have. The 
actual phase-space structure cannot be simpler than that. 

A third approach to studying the structure of galactic halos is by means of N-body sim- 
ulations Initial density perturbations are chosen according to a theoretically expected 
probability distribution, e.g. Gaussian density perturbations with the CDM power spectrum 
evolved from a flat Harrison-Zel'dovich primordial spectrum. Overdensities grow in the sim- 
ulations and form structures which are identified with galactic halos. In the large N limit, 
the outcome of the simulations can be identified with the predictions of CDM cosmology. 
However, present constraints on computing power restrict N to be of order 10 6 per galactic 
halo. In contrast, the number of axions per galactic halo is of order 10 84 , and the number of 
WIMPs of order 10 68 . Since phase-space is 6-dim. , the resolving power of the simulations is 
only a factor of ten in each dimension, including the three physical space dimensions. The 
simulations do show that galactic halos are not thermalized and that they contain discrete 
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flows. However, the limited resolution allows only a few flows to be seen. The simulations 
predict that a large galaxy like our own has hundreds of satellites. The tidal disruption of 



these satellites produces low velocity dispersion streams ]nj which are in effect localized 
flows. 

The purpose of this paper is to analyze the effect of the Sun's gravitational field on a 
flow of collisionless dark matter through the solar system. This effect was analyzed several 



years ago by K. Griest [2D| in the case of the isothermal model. He calculated the annual 
modulation of the signal in WIMP detectors including the effects both of the Earth's orbital 
motion and the Sun's gravity. Earlier discussions of the effect of the Sun's gravity on a flow 
of ordinary collisionful matter can be found in papers by H. Bondi and F. Hoyle |2I|| , J.M.A. 



Danby and G.L. Camm ||22|| , and Danby and T.A. Bray |[23|| . We will study mainly the effect 
of the Sun's gravity on a single flow of collisionless dark matter in the limit of zero velocity 
dispersion. The effect of finite velocity dispersion can be taken account by integrating over a 
distribution of initial velocities. Our results can be applied to the calculation of the annual 
modulation of the signals in both the axion and WIMP dark matter detectors. We plan 
to do this in detail in a subsequent paper which will be concerned with the experimental 
implications of the present investigation. In the remainder of this Introduction, we give a 
general discussion of the annual modulation of the signal in the axion and WIMP detectors, 
including the role of the Sun's gravity. 

A flow produces a peak in the energy spectrum of photons produced by axion — ► photon 
conversion in the cavity detector of dark matter axions 0. The peak's frequency uo is 
determined by the energy conservation relation, Tiu = m a (c 2 + |ffl ow ^i a b) > i n terms of the flow 
speed relative to the laboratory. Because the Earth rotates, there is a diurnal modulation of 
the peak's frequency. By measuring the average frequency and the amplitude and phase of 
the modulation, one can derive the velocity vector -Ufl ow ^_ e of the flow relative to the Earth 
(©). The height of the peak is proportional to the flow density. 

On top of the diurnal modulation, there is a much larger annual modulation because the 
Earth orbits the Sun (©): 

hu)(t) = m a [c 2 + i({T flow ^e(i)) 2 ] = m a [c 2 + -(iW*-©(f) - v®+ Q (t)) 2 ] , (1-2) 

where u is the peak frequency averaged to remove the diurnal modulation. v®+ Q (t) follows 
from the Earth's orbital motion. Eq. (|1.2|) then gives u(t) in terms of v^^q, the flow 
velocity relative to the Sun. However, galactic halo models give the velocity of the flow 
long before it reaches the Sun. For example, for the 'big flow' of Eq. ( [LTD , ^oo = v 5 — v Q 
where v Q is the velocity of the Sun relative to the Galaxy. ■Wa ow ^ differs from because 
the flow is affected by the gravitational field of the Sun. The relationship between and 
iJfiow^o depends on position and hence, as the Earth orbits the Sun, on time. 

One of the main goals of this paper is to derive this relationship. Note that a single flow 
of velocity far from the Sun, produces more than one flow in the vicinity of the Sun. We 
call n(f) the number of flows at position r. We want to calculate Vq ow ^q > j(f) for each of the 
flows j = l...n(f). We also want to calculate the density dj(f) of each of the flows. 

The simplest case is when the Sun is treated as a point mass. Then, n = 2 everywhere. 
Indeed there are two different ways in which particles coming with velocity Voo from far away 
can reach an arbitrary location r. See Fig. 1. We will find a caustic line downstream of 
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the Sun, on the positive z-axis if z is the unit vector in the direction of Voo. Indeed, every 
trajectory passes by some point on the positive z-axis; the latter is a focal line. We will 
see that the density diverges there as - where p is the distance to the z-axis. We call this 
caustic the 'spike'. 

In reality the Sun has finite radius. This implies that there is a maximum scattering 
angle GmaxC^oo) < 180°. As a result, at large distance r from the Sun, n = 1 for 6 > max 
and n = 3 for 6 < max , where 6 is the polar angle relative to the direction of Voo = zv^. On 
the conical surface which separates the n — 1 and n = 3 regions, there is a caustic, which 
we call the 'skirt'. The spike caustic along the z-axis downstream from the Sun remains 
as in the point mass case. As the Earth moves around the Sun, it may come close to the 
spikes associated with some of the flows and it will go through many of their skirts. When 
it approaches these caustics, the density becomes larger. We will see that the profile of the 
skirt caustic is a functional of the mass distribution inside the Sun. We may optimistically 
envisage a time where the axion dark matter detector on Earth, monitoring the peak sizes 
associated with various flows, provides information about the solar mass distribution. 

WIMP detectors measure the nuclear recoil energy when a WIMP scatters elastically off 
a nucleus in the detector |||J. The events caused by a single flow of WIMPs produce a 
plateau in the distribution of measured recoil energies. The edge of the plateau is at the 
highest recoil energy compatible with the kinetic energy of the WIMP. The highest recoil 
energy occurs when the collision is back-to-back in the center of mass. It is therefore a 
simple function of the mass of the nucleus, the mass of the WIMP and the speed of the 
WIMP in the laboratory frame. In principle, a WIMP detector can measure the speed of 
each of the flows by measuring the edge of the corresponding plateau in the recoil energy 
distribution. If the detector has directionality [|4j], it may determine the directions of the 
flows as well as their speeds. 

All WIMP detectors are presently sensitive to the local dark matter velocity distribution 
through the annual modulation effect || . The event rate in a WIMP detector is proportional 
to the flux of WIMPs through the detector. (In contrast, the axion to photon conversion 
rate in the cavity detector of dark matter axions is proportional to the density of axions, 
not their flux.) The WIMP flux changes in the course of the year as the velocity of the 
Earth's orbital motion around the Sun adds to or subtracts from the velocity of the Sun 
with respect to the galactic halo. 

Interestingly, the annual modulation effect in WIMP detectors distinguishes between the 
isothermal and caustic ring models of the Galaxy. In the isothermal model, the halo has 
zero average velocity everywhere. It has no net angular momentum, no rotation. The Sun 
moves with velocity of order 220 km/s in the direction <pG of the galactic disk rotation. The 
Earth moves with speed 30 km/s around the Sun in a plane, called the ecliptic, inclined at 
approximately 60° relative to (pa- The direction of is closest to (fie on about June 2. The 
WIMP flux, and hence the signal in a WIMP detector, is highest then. It is lowest six months 
later, near December 2. The amplitude of the modulation is of order cos(60°) 2 3 2 ° Q ~ 7%. 
In the caustic ring model, on the other hand, the halo does have net angular momentum. 
The locally dominant flows, such as Eq. (|1 . 1| ) , have velocity components in the <pc direction 
of order 470 km/s, more than twice the velocity of the Sun in that direction. Relative to 
the Sun, the average flow of dark matter is in the +4>g direction in the caustic ring model, 
whereas it is in the — 4>g direction in the isothermal model. So the annual modulation in the 
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caustic ring model |14|p5| -[27| has opposite sign to that of the isothermal model: the event 



rate in WIMP detectors is largest near December 2, lowest near June 2. Non-isothermal 



galactic halo models other than the caustic ring model are discussed in ref. |28 



The DAMA experiment in the Gran Sasso laboratory claims to have observed the annual 



modulation of the WIMP signal p9| . The DAMA results are consistent with the prediction 
of the isothermal model, and inconsistent with the caustic ring model. However, almost all 
of the region of WIMP mass and elastic scattering cross-section space implied by the DAMA 
signal is ruled out by the null result of the CDMS dark matter search |30|. Some of the 
region is also ruled out by the null result of the EDELWEISS experiment |5TJ]. 

At any rate, it is clear that the annual modulation is an important way to distinguish 
signal from background in WIMP searches, and that it deserves the best possible theoretical 
treatment. Most calculations and discussions neglect the effect of the Sun's gravity. However, 
as we mentioned already, the Sun's gravity effect was analyzed by K. Griest for the case of 
the WIMP signal in the isothermal model |[20|| . In that model the flow of dark matter relative 
to the Sun produces a wake in the — 0g direction. We may think of it as a spike caustic 
except that it is spread out in space because of the large velocity dispersion (ViP ~ 270 
km/s) of the flow. We reanalyze the effect of the Sun's gravity on the annual modulation of 
the WIMP signal in section V. Our results there are in qualitative agreement with those of 
Griest. 

In summary, our purpose is to study the effect of Sun's gravity on a flow of collisionless 
dark matter in the limit of zero velocity dispersion. The effect of finite velocity dispersion 
can be taken account by integrating over a distribution of initial velocities. The results can 
be applied to the flows, or streams, predicted by any specific halo model. 

This paper is organized as follows. In section II, we give a general discussion of cold 
(i.e., zero velocity dispersion) flows and associated caustics. We derive formulas for the 
density of such flows that will be useful in subsequent sections. In section III, we derive the 
density and velocity distributions of a stationary cold dark matter flow past a point mass. 
In that case, all the results are analytical and exact. In section IV, we derive the density 
and velocity distributions of a stationary cold dark matter flow past a spherically symmetric 
mass distribution like that of our Sun. In this case, the trajectories that go through the 
Sun are calculated numerically. The skirt caustic phenomenon appears here. In section V, 
we extend our formalism to dark matter velocities with arbitrary velocity distribution, and 
apply the results to the isothermal model as an example. In section VI, we summarize our 
findings. 

II. DENSITY FORMULAS 

Here we derive expressions for the physical space density of a flow with zero velocity 
dispersion, using the general approach of ref. ||12|| . In such a flow each particle carries a 3- 
dimensional label a = (a\, a 2 , a 3 ). The flow is entirely defined by giving the positions x(a, t) 
of all particles at all times. The velocity of particle a is v(a,t) = dx(a,t)/dt. In general, 
the map a — > x is one to many. Hence, at each point r of space, there is a discrete set of 
flows with velocities Vj(r, t) = v(dj(r, t),t), j = l...n(f, t), where <Sj{r, t) are the solutions of 
r — x(a,t). n(f,t) is the number of distinct flows at (r,t). The total number of particles is 
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N= d 3 a 



d 3 N 



da,\da,<2,doiz 



la 



J ' ^ T ^{da 1 



d 3 N 



da.ida.-i, 



det (1 



The density of particles in physical space is therefore: 



where 



dm = j2 



d 3 N 



' da\da 2 da^ 



{aj{f,t)) 



(Ox \ 
as)' 



ct=ctj (f,t) 



D(d,(f,t),t)\ ' 



(2.1) 



(2.2) 



(2.3) 



The formula for the density is manifestly a — *> (3(a) reparametrization invariant. 

The number of flows n(f, t) can change from point to point. Wherever a change occurs, n 
changes necessarily by two, i.e. two flows are added or two flows disappear. At those places 
the map a —>■ x is singular, D vanishes and the density d is infinite. Such places are called 
caustics. Caustics are generically surfaces because the condition D = defines a surface. 
Only in degenerate cases are caustics lines or points. 

A convenient generic parametrization of a zero velocity dispersion flow is as follows. 
Choose a surface S that each particle crosses at least once in the course of its trajectory. S 
may be time-dependent. Let a = (ai,a 2 ,t ) where t is the time the particle crosses S for 
the first time, and («i, a 2 ) labels the crossing point on <S. Then 



d 3 N 



d 3 N 



da 3 da\da 2 dt^ 



F(a 1 ,a 2 ,t ) 



d 2 S 
da\da 2 



(a ll a 2) t ) 



(2.4) 



where _F(ai,a2,£o) is the flux of particles through the surface. For a stationary flow it is 
natural to choose a time-independent surface S, and a time-independent parametrization 
(ai,«2) of the points on S. F and d ^ a2 are ^-independent then, and x(a\, a 2 , t ; t) = 
x(a 1 ,a 2 ,t - 1 ). 

In the case of a flow which is initially uniform in space, with density ^^(to) arid velocity 
Voo = Voo(to)z, and which is incident upon a spherically symmetric mass distribution, the 
flow is axially symmetric about z. It is natural to choose for S a plane perpendicular to 
z at large negative z and (01,0:2) = {b, ip) where b is the impact parameter and 95 is the 
azimuthal angle about the f-axis. In this case: 



x(b,if,t ;t) = (x cosip + y sin <p) p(b, t ; t) +zz(b,t ;t) 



(2.5) 



provided we adopt the convention that p can be negative as well as positive and restrict 
< ip < 7r. We also let b have either sign. [In the usual cylindrical coordinates, where p > 
and < ip < 2tc, Eq. ( |2.5|) would be invalid because ip changes by 7r when the z-axis is 
crossed.] From Eq. (|2.5|) follows: 




p det 



' d(p,z) 
d(b,t ] 



(2.6) 
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Moreover 



da 3 

Hence, Eq. fl2.2|) becomes 



■' :V FM* S 



dbdip 



doo(t ) i>oo(*o) |6| 



„ x = J_ y |&|<M*o)Uoo(*o) 



Idet 



(2.7) 



(2.8) 



b=bj (p,z,t), t =t 0j (p,z,t) 



where bj(p, z, t) and toj(p, z, t) are the solutions of p{b, t Q ; t) = p, z(b, to] t) = z. If the flow is 
stationary, d^ and are t ~i n dependent, p(b, t ; t) = p(b, t—t ), and z(b, t ; t) = z(b, t—t ). 
In that case 



d(p,z) 



d v n 



E 



\b\ 



det (wi) 



(2.9) 



b=bj (p.z). t=tj (p,z) 



where the tj and bj are the solutions of p(b, t — to) = p and z(b, t — t ) = z. 



III. COLD FLOW PAST A POINT MASS 

In this section we derive the density and velocity distribution of collisionless particles 
in a stationary zero velocity dispersion flow past a point mass M. The flow is axially 
symmetric as well as stationary. Let doo and Voo = v^z be the uniform density and velocity 
in the absence of M, and (z, p, ip) cylindrical coordinates centered on the point mass. The 
trajectory of the particle with impact parameter b is given in parametric form by 



t(b,ip) = — (e sinh ip — ip) 

z(b,tp) = a{e sinh^ + 1 — - expip) 

p(b,i/>) = b(l - -expV), (3.1) 
e 



where a = ip is the "eccentric anomaly" parameter, and e = y 1 + (^j . Fig. || illustrates 
a pair of trajectories. As in the previous section, we do not use cylindrical coordinates in 
the conventional way. We allow p to have either sign, but restrict ip to take values between 
and 7r. Also we let b have either sign. 

The velocity components of the flow are: 



v p = — — = -v , 



dp dip b exp ip 

•rn ib 

(3.2) 



dtp dt ea(e coship — 1) 

dz dip e cosh^ — - expip 



dip dt " e cosh^ — 1 



The density distribution is given by Eq. (|2.9|) . In the point mass case, n = 2 everywhere. 
Upon changing variables t — > ip, we obtain: 



9 



7, ^OO ^OO \ ^ 

3=1 



\P\ 



\b\ 


dt 
dip 




det ( 


d(p,z) \ 
d(b^)) 



b=b j (p,z),4!=ip J (p,z) 

(ecoship — 1) 



\ (1 - iexpV0p2(&,V0l 



(3.3) 



b=b j (p,z),ip=ip j (p,z) 



where bj(p, z) and if)j(p, z) are the two solutions of p(b, if)) = p and z(b, if)) = z, and 

= ^ ( e exp(— *0) - 1 + (e — expip + exp(2^)^ . 
Eqs. ( |3.1| ) can be inverted: 

b±(p,z) = £ 



(3.4) 



1 ± Jl + V 



e±~ ( I T v' l '/ 



(3.5) 



where e± = y 1 + (^) 2 and 



4a . . 
y = — (r + z) 

P 



2a 



r sm 2 (6/2) 



(3.6) 



with r = y/ p 2 + z' 2 and = tan l (p/z). Henceforth we label the two flows with j = ±. 
Inserting Eqs. ( |3.5|) into Eqs. (|3.4j) and ( |3~3"1) yields: 



£ 2 ± = r 



(3.7) 



and 



± 2 



Note that d + — gL = doo everywhere. The velocities are: 



w P ± = - — (6± - p) 

r 



V t ± = Vgo [l + 



a p 
r b+ 



The total density 

d = d+ + d_ 



d, 



oo 

2~ 



2a 



r sin 2 (#/2) 



1 + 



2a 



r sin 2 (0/2) 



(3.8) 



(3.9) 



(3.10) 
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It diverges on the positive z-axis as 



V2az 

a ^ — — — . (3.11) 

H 

This divergence occurs because the flow is focussed onto the positive i-axis by the point 
mass' gravity. Each trajectory crosses the positive 5-axis, at z = Jr. The positive i-axis is 



the location of a line caustic, which we call the 'spike'. Eq. ( |3.11| ) is valid also for a finite 



size but spherically symmetric mass, at z large enough that the trajectories there did not 
pass through the mass. 

As we remarked already in section II, caustics are generally surfaces rather than lines 
or points A line or point caustic is a degenerate case. The caustic of Eq. ( J3.11| ) is a 



line only because of the axial symmetry of the flow. If less symmetry is present - e.g. if the 
mass has a quadrupole moment, or if the initial velocity field of the flow is inhomogeneous 
- there is still a caustic downstream of the mass but it is a surface. Since it can collapse to 
a line, the surface is a tube of some sort. The structure of the tube is a problem left for 
future investigations. 



IV. COLD FLOW PAST A SPHERICALLY SYMMETRIC MASS 

In this section we derive the density and velocity distribution of a zero velocity dispersion 
flow past a spherically symmetric mass distribution of total mass M®. The mass distribution 
is characterized by a radius R® outside of which the density rapidly approaches zero. The 
trajectories passing through the mass must be calculated numerically. As before, the only 
force acting is the gravity of M®. The density and velocity = v^z of the flow long 
before encountering the mass distribution are assumed to be constant in both space and 
time. Axial symmetry is preserved and the flow is stationary. 

The density and velocity distributions derived in the previous section ( |111| ) are valid 
here provided the trajectories did not pass through the mass distribution at any point in 
their past. Defining b® to be the impact parameter of the trajectory grazing the surface 
of the mass distribution, all results of the point mass case hold for trajectories labeled by 
b > b®. The b > b® trajectories are given in parametric form by Eqs. (|3.1|) which implies 
r = a(e cosh ip — 1). The point of closest approach to the mass is where ip = 0. Setting ip = 
and r = R® gives the grazing impact parameter 

b® = R®^l + 2a/R®. (4.1) 

On the other hand, all trajectories for b < b® will pass through the mass distribution at 
some point in their motion. 

Here it is convenient to work in spherical coordinates (r, 9, <p). As in the previous section 
the azimuth (p is restricted to be between and it. We measure the polar angle 9 from the 
+z-direction and it takes values between —n and +7r. In this system of coordinates the 
polar angle smoothly changes sign as a trajectory crosses the symmetry axis downstream 
from the mass. As before, we let b have either sign. 

We derive the density distribution of a stationary flow with axial symmetry in spherical 
coordinates. Analogously to Eqs. fl2.6| ) and Q2.9p, we have 
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- 2 sin 9 det 



' d(r, 9) ' 



(4.2) 



and 



d(r,9) 



r 2 | sin#| 



E 



|6| 



|det 



(4.3) 



b=bj(r,6),t=tj(r,t 



The index j counts the flows. In the point mass case, particles with impact parameter close 
to zero are scattered by an angle close to 180°. In that case, n = 2 everywhere. On the 
other hand, if the mass has finite size, there is a maximum scattering angle |0 m ax| < 180°. 
This angle defines a cone. Far from the mass, n = 1 for \9\ > |0 max | whereas n = 3 for 
\9\ < |0 max |. Since the number of flows changes on the surface of the cone, we expect a 
caustic there. 

We break the sum over j into two parts. The first sum is over the flows for which we 
have calculated an exact result: the b > b Q flows at any time in their motion and the b < b & 
flows before entering the mass distribution. We label these flows by the index k (= 1, n&). 
The second sum is over the flows which must be calculated numerically, namely the b <b Q 
flows at any time after they have entered the mass distribution. We label these flows by the 
index I (= 1, ...,n{). The density is then 



hi 



d(r,f) = E*M) + #iS 

r 2 |sm#| fri 



\b\ 



I det fM) 
| aet \d(b,t)J 



(4.4) 



b=b l (r,9),t=ti{r,( 



where the dk{r, 9) expressions were derived in section |Tl| and the terms in the sum over I we 
address next. Although we have not shown it explicitly Eq. ( |4.4j ), n& and n\ are dependent 
upon (r, 6>). We will see that (n k ,rii) = (1, 0), (1, 2), (2, 1), (0, 1) or (0,3), depending on 
position. 

To simplify the determinant, while exposing the caustic features intrinsic to Eq. (14), 
we take the limit where r >> R & , a. In that case, for the Z-flows 



det 



'd(r,6y 



d9_ 

'db 



(4.5) 



Hence 



d,(r, 



doo 


M 




sin 9 




de 



(4.6) 



and the numerical computation of the second term in Eq. ( f4.4| ) becomes straightforward. 
For the /-flows at r >> R & and a, \9{b)\ is the scattering angle for the trajectory of impact 
parameter The scattering angle of a particle incident at impact parameter b with initial 
kinetic energy E upon a central potential V(u) is |32| 







7T 



|6| du 

y/l - V{u)/E 



b 2 u 2 



(4.7) 



12 



where u = l/r and u max = l/r min . We express the potential as V(x) = GMq<P(x)/Rq with 
x = r/R & . We choose the potential of the mass distribution to be 



[x = - 



a + (3x [ 



6 



1 + 7X 2 + fix 7 



(4.8) 



setting a = 5, (3 = 5000, and 7 = 8 to give a close fit to the solar mass distribution 
Thus, all of the following calculations should be considered relevant to the Sun. In Eq. ( 
~ —1/x for x > 1, the point mass limit, and ~ —a + 7X 2 for x << 1 reflecting the 
fact that the density in the solar core is constant. The potential ( f4.8| ) is plotted in Fig. ^| in 
comparison to = —1/x. 

We calculate the scattering angle 0(|6|) for the potential of Eq. ( ^4.8| ). For comparison, 
the scattering angle for the = —1/x potential is the exact result 

6(6) = - 2 arcsin (-^j sign(6) , (4.9) 

in which 0(6) can take values between —it and +tt. We plot 0(6) for the distributed mass 
and point mass cases in Fig. || calculated for two flow velocities, = 200 km/sec and 500 
km/sec. 

In Fig. H there is a maximum scattering angle max < tt- For < 9 < max , there 
are two flows labeled by different values of 6 with sign(6) = — sign(#). In addition, there is a 
flow for which 6 and 9 have the same sign. For large r, the third flow is in the sum over k. In 
the notation of section pTT|, it has impact parameter 6 + . So n = 3, for large r and |6*| < max . 
For max < \9\ < n, the two flows of Fig. |3| are absent but the 6 + flow is still present, so 
n = 1 there. 

Near the cone |#| = max is a caustic surface which we call the "skirt." The determinant 
in Eq. ( |4.4|) vanishes there since d6/db = 0. We calculated the maximum scattering angle 
0ma X as a function of flow velocity Voo ranging from 100 km/sec to 1000 km/sec. See Fig. [| 
for the result. 

We now examine the density Eq. ( |4.4| ) near the skirt caustic for angles 9 < max . The 
function 0(6) can be approximated near max as 

0(6)~0 max -i#(6-6 c ) 2 (4.10) 

where 6 C is the impact parameter for the maximum scattering angle, and H = — ^-\ b=bc . 



For 9 near to but smaller than max , we have = 1 and n\ = 2 in Eq. ( [4.4| ). The k = 1 
flow is the one labeled + in section [HI]. Using Eqs. ( |4.5| ) and ( |4.10| ), we get the approximate 
density near the skirt 



4 kirt ^ d + + ^tt^^ H ( eL-g) for ^ <0 -x (4.11) 
= d + for #>0 max . (4.12) 

A generic feature of caustic surfaces is evident in Eq. ( |4. 12| ) , namely that the density diverges 
near the surface as l/V^max; — 9 when the skirt is approached from 9 < max side. The 
density remains constant (equal to d + ) when the surface is approached from 9 > max . 
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The strength of the caustic is proportional to b c /r 2 sin max vi?. We plot H and b c over a 
range of flow velocities in Figs. [| and || respectively. For comparison, we plot b c alongside 
the grazing impact parameter 6 . 

Moving away from the skirt toward smaller values of \9\, an angle will be reached 
where the flows change from rik — 1, ni — 2 to rik = 2, ni = 1. 6 = 0(6©) is the scattering 
angle of the grazing trajectories. Using Eq. (|4.1|) gives 



6 = - 2 arcsin — — sign(6 ) . (4.13) 
\a + RqJ 

For |#| < 1 00 1, only one flow passed through the solar interior. We plot in Fig. | as a 
function of the flow velocity Voo. 



For r >> i? and a 



'max 



(n k ,m) = (1,0) for \9\ > |0 E 

= (1,2) for |0 | < \9\ < |0 max | (4.14) 
= (2,1) for \9\ < |0 | . 

There are two regions in or near the Sun where different n{) values arise. First, n{) = 
(0, 1) on the interior of the Sun on the upstream side. Second, there is a region near the 
Sun on the downstream side where (rik, n i) — (0, 3). The flow properties in these regions are 
calculable with Eq. (|4.4| ). We leave that task for future investigations. 

The velocity distributions for the fc-flows are given in Eq. ( |3.9| ). For the /-flows, using 
energy and angular momentum conservation, 



2a (bi\ 2 



r 



v 0l = Voo — . (4.15) 
r 

For the /-flows at large r, 9 is the scattering angle 0. Hence, the bi(9) are given by Fig. [|. 

Using 0(6) it is straight-forward to calculate the full expression for the density distri- 
bution at large r including the flows through the Sun. This result is given in Fig. [7] for the 
mean Earth-Sun distance. The density field of a b < b & flow, at large r, is smaller than the 

b 2 

sum of the b > b & flows, roughly by the factor -f . Away from the caustic skirt, the flows 
through the Sun contribute only 4 x 10 -3 to 4 x 10~ 5 of the total density at the Earth's 
location for velocities ranging from = 100 km/sec to 1000 km/sec. 



The limit Rq — ► seems paradoxical. In section |T| we found that n = 2 everywhere, 
whereas here n — 1 for |0| > |0 max | and n = 3 for \9\ < |0 max |- What happens is that when 
Rq — > and max — > n one of the three flows acquires vanishing density. The flow that 
disappears is the one with the smallest impact parameter for given scattering angle in Fig. || 
The flows that survive are the flow with the largest impact parameter for given scattering 
angle in Fig. |3| and the flow with sign(#)=sign(6), i.e., the flow labeled + in section |T|. 

To conclude this section, we briefly discuss the relationship between the density distri- 
bution formulas for a flow past a point mass and the Rutherford differential cross-section 
formula. The differential scattering cross-section is 



14 



da 



(0) 



sin 9 



db 
d6 



(4.16) 



and the Rutherford formula for scattering in a r 1 potential is 

dQ { )R 4sin 4 (0/2)' 



(4.17) 



where a = GM/v 2 ^ as before. Consider only the "— flow," in the notation of section |T|, 
in the limit where r » a and 9 ^ 0. As shown above, the determinant then simplifies, 

and the density can be written in terms of the Rutherford formula 



uel V d(b,t) 

3q. O: 



d_(r,0) 



dn 



r 2 sin 6 1 



de 



d^ da 

~^dn 



(0)r 



(4.18) 



Using a fact derived in section [TTT| , namely that d + — oL = cf^ everywhere, we get 



d+(r,0) =41 + 



4r 2 sin 4 (0/2) 



(4.19) 



One can verify that the Eqs. ( f4.18p and ( |4.19| ) agree with the exact results for d± derived 
in section |T|. Taking the r — ► oo limit of Eq. ( |3.8|) at fixed 9 ^ and keeping the terms to 
leading order yields Eqs. (4.18|) and (|4.19). 



V. HOT FLOW PAST A POINT MASS 

In many discussions of dark matter detection on Earth, it is assumed that the distribution 
in the galactic halo is isothermal. The local velocity distribution is then 



>7ra r cr 



in the absence of the Sun, for an observer moving with the Sun's velocity v & . a = 
y^l < v 2 > ~ 220 km/s, where < v 2 > x ^ 2 is the velocity dispersion of the galactic halo. 
a and v® are approximately equal. As was mentioned in the Introduction, there is little 
theoretical motivation for assuming an isothermal distribution. Indeed most of our paper 
is concerned with discrete flows, which we believe to be the relevant description. However 
our methods can also be applied to smooth velocity distributions such as Eq. ( |5.1|) . This is 
done in this section. We neglect here the Sun's size, treating it as a point mass. 

If the velocity distribution in the absence of the Sun is f(voo) independently of position, 
what is the position-dependent velocity distribution f Q (r,v) when the effect of the Sun's 
gravity is included? The answer to this seemingly daunting question is p2]j20 



f Q (f,v) = f(v 00 (r,v)) (5.2) 

where ■y 00 (^*, v) is the initial velocity of the trajectory that passes by position r*with velocity 
v. Indeed, Liouville's theorem states that phase-space density is conserved along particle 
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trajectories. Since f(voo) is the initial phase-space density, when the particles are far from 
the Sun, and f Q (r,v) is the final phase-space density, Eq. ( |5.2|) follows. 
To obtain v^r, v), we exploit the constants of the motion: 

2 r 



I = f x v 

A = v x (r x v) — GM Q f 



(5.3) 



E, I, and A are the energy, angular momentum and Laplace-Runge-Lenz vector, per unit 
mass in each case. In the initial state, 



E= ~ 2 vL 

A = bv^p + GM Q z. 



(5.4) 



where b is the impact parameter, z is in the direction of the initial velocity (voc = Vooz), ip is 
perpendicular to the plane of the trajectory and p = (p x z. The notation is consistent with 
that of section III. Since 



2GM, 







(5.5) 



it only remains to express z in terms of r and v. We have: 

A x 1 = +bvl D (ap - bz) 
where a = ^ir 2 -, as before. Hence 



Substituting Eqs. ( |5.3[ ) into Eq. ( {5.7Q , we obtain after some algebra: 



1 



tV, r,t> 



I 2 — av 2 r — av^ir ■ v) 



+ rav„ 



1^ ^ v 1 
-(r ■ v) H 

r r 



a?vl Q + I 2 

where I 2 = r 2 v 2 — (r • v) 2 and is given in terms of r and v by Eq. (|5.5| ). 



If the isothermal velocity distribution (5.1) is assumed, we have therefore: 



fe(r,v) 



(In 



e M--t v G(r, v)]©(t&(r, *>)) 



where 0(x) = for x < 0, 1 for x > 0, and 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



«g(^ *0 = + ^oo(r, w)) 2 = v% + + 2w 



v cos a 



I 2 — av^r — avoo(f- v) + av 2 ^ cos (3 



2 

_ _, v r 

r ■ v H av c 



(5.10) 
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where a is the angle between v e and the dark matter velocity v and (3 is the angle between 
v & and the position r of the observer relative to the Sun. A formula for v q(t, v) was obtained 



by J.M.A. Danby and G.L. Camm |22| over 40 years ago. K. Griest |20[ obtained a formula 
which differs from that of Danby and Camm. Our formula appears different from both that 
of Danby and Camm, and that of Griest. We verified that the RHS of Eq. ( |5.8| ) is indeed 
v 2 ^, and that Eq. (|5.8|) checks out in the I = and / — > oo limits, in the a = limit, and 
when Voo \\ r. 

The annual modulation of the event rate in WIMP dark matter searches is caused for the 
most part by the seasonal variation of the flux of dark matter particles as the Earth's orbital 
motion adds to or subtracts from the motion of the Sun through the galactic halo. As was 
discussed in the Introduction, there is an additional contribution because the Sun's gravity 
modifies the density and velocity distributions in the solar neighborhood. The combined 
effects of the Earth's orbital motion and the Sun's gravity can be computed if we know 
the distribution of dark matter in the Earth's frame of reference. Following Griest, that is 
obtained by simply replacing v Q in Eqs. ( |5.9|) and ( |5.10| ) by v Q + v e where % is the velocity 
of the Earth relative to the Sun. The event rate depends also on the WIMP mass and 
such experimental variables as the detector threshold, the target mass and, of course, the 
scattering cross-section. 

The effect of the Sun's gravity on the event rate in WIMP searches is largely due to the 
change in the density of dark matter. Moreover the signal from axion to photon conversion in 
the cavity detector of galactic halo axions is simply proportional to the dark matter density. 
So let us calculate the density d(r, 9) where r is the distance to the Sun and 9 is the angle 



between r and v Q . We could integrate Eq. ( |5.9| ) over v. However, it is easier to integrate 

d(r,9) = [ dV/^ooM^iU (5.11) 



where d(r, v^) is the density due to a single flow of initial velocity v^, Eq. ( |3.10|) . This 
yields 



^ poo /■2-7T 

d(r,9) = / u 2 du I ds dip exp[— (u 2 + u 2 G) + 2uu Q s)} 

7T2 JO J-l JO 




(5.12) 



where m r = 



4 



X — — — — m 2 (1 — s cos 6 — y/l — s 2 sin 6 cos (p) (5.13) 



and a a = ^sa. The integral in Eq. fl5T2D was evaluated numerically. Fig. || shows d(r) in 



case the Sun were at rest with respect to the halo (v Q = 0). At the Earth's distance from 
the Sun, the density increase due to the Sun's gravity is approximately 3%. 

Fig. ^| shows d(r, 9) as a function of 9 for v & = a = 220 km/s and r = 1.5 x 10 13 cm, 
the mean Earth-Sun distance. At the Earth's distance from the Sun, the density increase 
is approximately 0.5% upstream {9 = n) and 6% downstream {9 = 0). The large density 
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downstream is the spike caustic but spread out because of the large velocity dispersion. The 
Earth's orbit is inclined at approximately 60° relative to v Q , with the Earth being most 
closely downstream of the Sun on March 3. Fig. |9] shows that the density variation between 
September 2 (8 = 120°) and March 3 {9 = 60°) is of order 1.4%. This is consistent with the 
earlier calculation of Griest p0[ . 



VI. CONCLUSIONS 

We analyzed the effect of the Sun's gravity on a flow of collisionless dark matter through 
the solar neighborhood. We were mainly concerned with the case of a cold flow, i.e. a flow 
with zero velocity dispersion. The cold flow results can be generalized to arbitrary flows by 
integrating over the initial velocity distribution. We neglected the effect of the Earth's 
gravity, but this is indeed small. The ratio of the Sun and Earth gravity effects is of order 
Tp-Tj- — 0.07 where M e is the Earth's mass, i? e its radius, and r is the mean Earth-Sun 
distance. 

In section II, we obtained formulas for the density and velocity distributions of cold flows 
with axial symmetry. In section III, these formulas were applied to a cold uniform flow past 
a point mass. In that case the number of flows n = 2 everywhere. There is a caustic line, 
which we called the 'spike', downstream of the Sun. It is located on the positive f-axis 
where z is the direction of the incoming flow (v^ = VooZ) . The profile of the spike is given in 
Eq. ( |3.11| ). The densities of the two flows everywhere are given by Eqs. ( |3.8| ), the velocity 



distributions by Eqs. (|3.9| ). These are exact results. 

In section III, we took account of the Sun's finite radius. We assumed that the solar 
mass distribution is spherically symmetric. Many of the results derived in section III are 
still valid. In particular, the spike caustic remains and its profile is still given by Eq. ( |3.11| ) 



for z large enough that the trajectories there did not go through the Sun in the past. A 
new caustic appears which we called the 'skirt'. It is located near the cone of axis z and 
opening angle O max where max is the maximum scattering angle. Inside this cone n = 3, 
whereas n = 1 outside the cone. The density enhancement associated with the skirt falls off 
as where r is the distance to the Sun, and is rather feeble on Earth even at the moment 
of skirt crossing. However, the profile of the density near the skirt depends on the mass 
distribution inside the Sun. In principle, by monitoring the densities of the various dark 
matter flows as we orbit the Sun, the mass distribution of the Sun can be investigated. The 
density profile near a skirt was obtained in terms of quantities H, b c and O max , which we 
defined and computed as a function of v oo for a realistic model of the solar mass distribution. 

In section IV, we applied our methods to the single flow of dark matter predicted by the 
isothermal model. That flow has direction — 0g opposite to the Sun's motion in the Galaxy 
and velocity dispersion — 270 km/s. We derived a formula which gives the velocity 
distribution at any point in the neighborhood of the Sun in terms of an arbitrary initial 
velocity distribution /(tw,). We also computed the density as a function of position relative 
to the Sun. There is a spike caustic downstream of the Sun but it is spread out because of 
the large velocity dispersion. In the isothermal model the density is largest (smallest) near 
March 3 (September 2). The amplitude of this density modulation is of order 0.7%. 
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FIG. 1. Two trajectories past a spherically symmetric mass distribution. The incident flow is 
assumed to be uniform and to have vanishing velocity dispersion. If the scatterer is a point mass, 
the particles with small impact parameter are scattered by an angle close to 180°, and hence the 
number of flows n = 2 everywhere. If the scatterer is distributed over a finite radius, there is a 
maximum scattering angle © max < 180°. In that case, n = 1 upstream of the scatterer, and n = 3 
downstream. 



22 








0.5 1 1.5 2 

x=r/R Q 

FIG. 2. A plot of the dimensionless Solar potential <fi(x) given in Eq. ([|]^) and compared with 
the point mass potential <p(x) = —1/x. 
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FIG. 3. The scattering angle for dark matter incident upon the Sun as a function of impact 
parameter, for flows of initial velocity 200 km/sec and 500 km/sec. The solid curves are for the 
distributed mass described in section IV and the dashed curves are for the point mass. The flows 
with lower initial velocities have larger scattering angle at large impact parameter. 
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FIG. 4. The maximum scattering angle @ max (solid line) and the scattering angle 0© of the 



grazing trajectory (dashed line) as a function of initial flow velocity v c 
scribed in section IV. 



for the solar model de- 
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FIG. 5. The parameter H = —^-\ b=bc as a function of the initial flow velocity t>oo, for 
solar model described in section IV. The strength of the skirt caustic is proportional to H~ 1 / 2 
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FIG. 6. The impact parameter b c of the trajectory of maximum scattering angle (solid line) and 
the impact parameter 6 of the grazing trajectory (dashed line) as a function of initial flow velocity 
v^, for the solar model described in section IV. The strength of the skirt caustic is proportional 
to b c . 
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FIG. 7. The density of a cold collisionless flow past the Sun as a function of polar angle for 
= 1 AU and Voo = 255 km/s. The divergence at 9 = is the spike caustic. The skirt caustic is 
6 = 115°. It is invisible on the scale of the main plot but is shown in the inset. 
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FIG. 8. The density as a function of distance to the Sun, in the isothermal model with 



if2. — 



270 km/s, if the Sun were not moving with respect to the halo (vq = 0). The increase in density 
at small r, compared to doo, is the effect of the Sun's gravity. At the Earth's distance the increase 
is approximately 3%. 
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FIG. 9. The density as a function of polar angle 9 at r = 1 A.U. in the isothermal model with 
ViP = 270 km/s and vq = 220 km/s. 9 is relative to the direction of —Vq. The Earth's orbit is 
inclined at approximately 60° relative to vq. The annual modulation due to the Sun's gravity can 
be read off by comparing 9 ~ 60° (March 3) and 9 ~ 120° (September 2). It is approximately 
0.7% in amplitude. 
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